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ABSTRACT 

We present the spectroscopic data for the galaxies studied in the EFAR project, which 
is designed to measure the properties and peculiar motions of early-type galaxies in 
two distant regions. We have obtained 1319 spectra of 714 early-type galaxies over 
(SJ ■ 33 observing runs on 10 different telescopes. We describe the observations and data 

reductions used to measure redshifts, velocity dispersions and the Mg6 and Mg 2 Lick 
linestrength indices. Detailed simulations and intercomparison of the large number 
of repeat observations lead to reliable error estimates for all quantities. The mea- 
surements from different observing runs are calibrated to a common zeropoint or 
^ e— | scale before being combined, yielding a total of 706 redshifts, 676 velocity dispersions, 

O 676 Mg6 linestrengths and 582 Mg2 linestrengths. The median estimated errors in 

the combined measurements are Acz=20kms _1 , Aa/a=9.1%, AMg6/Mg6=7.2% and 
AMg2=0.015 mag. Comparison of our measurements with published datasets shows 
no systematic errors in the redshifts or velocity dispersions and only small zeropoint 
corrections to bring our linestrengths onto the standard Lick system. We have assigned 
galaxies to physical clusters by examining the line-of-sight velocity distributions based 
on EFAR and ZCAT redshifts, together with the projected distributions on the sky. We 
derive mean redshifts and velocity dispersions for these clusters, which will be used 
■ in estimating distances and peculiar velocities and to test for trends in the galaxy 

population with cluster mass. The spectroscopic parameters presented here for 706 
galaxies combine high quality data, uniform reduction and measurement procedures, 
and detailed error analysis. They form the largest single set of velocity dispersions and 
linestrengths for early-type galaxies published to date. 

Key words: galaxies: clustering — galaxies: distances and redshifts — galaxies: 
elliptical and lenticular, cD — large scale structure of universe — surveys 
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1 INTRODUCTION 

We are measuring the peculiar motions of galaxy clus- 
ters in the Hercules-Corona Borealis (HCB) and Perseus- 
Pisces-Cetus (PPC) regions at distances between 6000 and 
15000 kms -1 using the global properties of elliptical galax- 
ies. This study (the EFAR project) has as primary goals: 

(i) characterising the intrinsic properties of elliptical galax- 
ies in clusters by compiling a large and homogeneous sam- 
ple with high-quality photometric and spectroscopic data; 

(ii) testing possible systematic errors, such as environmen- 



tal dependence, in existing elliptical galaxy distance estima- 
tors; (iii) deriving improved distance estimators based on 
a more comprehensive understanding of the properties of 
ellipticals and how these are affected by the cluster envi- 
ronment; and (iv) determining the peculiar velocity field in 
regions that are dynamically independent of the mass dis- 
tribution within 5000 kms -1 of our Galaxy in order to test 
whether the large-amplitude coherent flows seen locally are 
typical of bulk motions in the universe. 

The background and motivation of this work are dis- 
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cussed in Paper I of this series (Wegner et al. 1996), which 
also describes in detail the choice of regions to study, the 
sample of clusters and groups, and the selection procedure 
and selection functions of the programme galaxies. In ear- 
lier papers we reported the photoelectric photometry for 352 
programme galaxies which underpins the transformation of 
our CCD data to the standard R magnitude system (Colless 
et al. 1993), and described our technique for correcting for 
the effects of seeing on our estimates of length scales and 
surface brightnesses (Saglia et al. 1993). This paper (Pa- 
per II) describes the spectroscopic observations and gives 
redshifts, velocity dispersions and linestrength indices for 
the programme galaxies. The CCD imaging observations of 
these galaxies, and their photometric parameters, are de- 
scribed in Paper III (Saglia et al. 1997), while descriptions 
of the profile fitting techniques used to determine these pa- 
rameters (along with detailed simulations establishing the 
uncertainties and characterising the systematic errors) are 
given in Paper IV (Saglia et al. 1997). The Mg-cr relation 
and its implications are discussed in Paper V (Colless et al. 
1998). Subsequent papers in the series will explore other in- 
trinsic properties of the galaxies and their dependence on 
environment, derive an optimal distance estimator, and dis- 
cuss the peculiar motions of the clusters in each of our survey 
regions and their significance for models of the large-scale 
structure of the universe. 

The structure of the present paper is as follows. In §^ 
we describe the observations and reductions used in ob- 
taining the 1319 spectra in our dataset (1250 spectra for 
666 programme galaxies and 69 spectra for 48 calibration 
galaxies) and discuss the quality of the data. We explain 
the techniques by which redshifts, velocity dispersions and 
linestrength indices were estimated from the spectra in §|§| 
including the various corrections applied to the raw values. 
In §^ we describe the method used to combine data from 
different runs and evaluate the internal precision of our re- 
sults using the large number of repeat measurements in our 
dataset. We then give the final values of the spectroscopic 
parameters for each galaxy in our sample: we have redshifts 
for 706 galaxies, dispersions and Mgfe linestrengths for 676 
galaxies and Mg2 linestrengths for 582 galaxies. We compare 
our results to previous studies in the literature to obtain ex- 
ternal estimates of our random and systematic errors. In §^ 
we combine our redshifts with those from ZCAT in order to 
assign sample galaxies to physical clusters, and to estimate 
the mean redshifts and velocity dispersions of these clusters. 
Our conclusions are summarised in §|(| 

This paper presents the largest and most homogeneous 
sample of velocity dispersions and linestrengths for elliptical 
galaxies ever obtained. The precision of our measurements is 
sufficiently good to achieve the goal of measuring distances 
via the Fundamental Plane out to 15000 kms -1 . 



2 OBSERVATIONS 

The spectroscopic observations for the EFAR project were 
obtained over a period of seven years from 1986 to 1993 in a 
total of 33 observing runs on 10 different telescopes. In this 
section we describe the spectroscopic setups, the observing 
procedures, the quality of the spectra and the data-reduction 



techniques. Further detail on these points is given by Bagg- 
ley (1996). 



2.1 Spectroscopic Setups 

Table |l] gives the spectroscopic setup for each run, including 
the run number, date, telescope, spectrograph and detector, 
wavelength range, spectral dispersion (in A/pixel), effective 
resolution (in kms -1 ), and the effective aperture size. Note 
that two runs (116 and 130) produced no useful data and 
are included in Table |l| only for completeness. Three runs 
utilised fibre spectrographs: runs 127 and 133 used Argus on 
the CTIO 4m and run 131 used MEFOS on the ESO 3.6m. 
All the other runs employed longslit spectrographs, mostly 
on 2m-class telescopes (MDM Hiltner 2.4m, Isaac Newton 
2.5m, Kitt Peak 2.1m, Siding Spring 2.3m, Calar Alto 2.2m) 
although some 4m-class telescopes were also used (Kitt Peak 
4m, William Herschel 4m, the MMT) . 

The spectra from almost all runs span at least the wave- 
length range 5126-5603A, encompassing the Mglfe 5 174 A 
band and the Fel 5207A and 5269A features in the rest- 
frame for galaxies over the sample redshift range cz~6000- 
15000 kms -1 . The exceptions are the spectra from runs 115 
and 131. Run 115 comprises 8 spectra obtained at the WHT 
with the blue channel of the ISIS spectrograph which have 
a red wavelength limit of 4970A (i.e. including H/3 but not 
Mgfo) . Since we have spectra for all these galaxies from other 
runs we do not use the redshifts and dispersions from run 
115. Run 131 comprises 128 spectra obtained at the ESO 
3.6m with the MEFOS fibre spectrograph to a red limit of 
5468A, including Mgb and Fel 5207A over the redshift range 
of interest, but not Fel 5269A beyond cz^llOOO kms -1 . For 
most of the runs the spectra also encompass H/3, and several 
span the whole range from Cal H+K 3933+3969A to Nal D 
5892A. 

The effective instrumental resolution of the spectra, Oi, 
was measure d fr om the autocorrelation of stellar template 
spectra (see §3.1 below), and ranged from 80 to 170 kms -1 , 
with a median value of 125 kms -1 . Both longslit and circu- 



lar entrance apertures were used. Slits were typically 1.7- 
2.0 arcsec wide and the spectra were extracted to the point 
where the galaxy fell to about 10% of its peak value. Cir- 
cular apertures (in the fibre spectrographs and the MMT 
Big Blue spectrograph) were between 1.9 and 2.6 arcsec in 
diameter. Further details of the observing setup for each 
telescope/instrument combination are given in Appendix A. 



2.2 Observing Procedures 

The total integration times on programme galaxies var- 
ied considerably depending on telescope aperture, observ- 
ing conditions and the magnitude and surface brightness of 
the target (our programme galaxies have R band total mag- 
nitudes in the range 10-16). On 2m-class telescopes (with 
which the bulk of the spectroscopy was done), exposure 
times were usually in the range 30-60 min, with a median of 
40 min; on 4m-class telescopes, exposure times were gener- 
ally 15-20 min (up to 60 min for the faintest galaxies) with 
single-object slit spectrographs, but 60 or 120 min with the 
fibre spectrographs (where the aim was high S/N and com- 
pleteness). Slit orientations were not generally aligned with 
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Table 1. The Spectroscopic Observing Runs 



Run 


Date 


Tele- a 


Spectrograph 


Detector 


A A 


AA b 


<Ji c 


Aperture d 




# 




scope 


+ grating 




(A) 


(A/pix) 


(km/s) 


(arcsec) 




101 


86/12 


MG24 


MarkIIIa+600B 


GEC 


4912- 


-6219 


2.27 


145 


1.9x 10% 


22 


102 


87/03 


MG24 


MarkIIIa+600B 


RCA 


4787- 


-6360 


3.10 


115 


1.9x 10% 


58 


103 


87/05 


MG24 


MarkIIIa+600B 


RCA 


4809- 


-6364 


3.07 


145 


1 Qv 10% 


37 


104 


88/04 


MG24 


MarkIIIa+600V 


RCA 


5025- 


-6500 


2.90 


125 


1.9x10% 


12 


105 


88/06 


MG24 


MarkIIIa+600V 


RCA 


5055- 


-6529 


2.90 


125 


1.9x10% 


37 


106 


88/09 


MG24 


MarkIIIa+600V 


Thompson 


5041- 


-6303 


2.21 


130 


1.9x 10% 


23 


107 


88/10 


MG24 


MarkIIIa+600V 


RCA 


5048- 


-6522 


2.90 


130 


1.9x 10% 


27 


108 


88/07 


MMTB 


BigBlue+300B 


Reticon 


3700- 


-7200 


1.14 


135 


2.5 


10 


109 


88/11 


KP4M 


RC+UV-Fast+17B 


TI2 


4890- 


-5738 


1.07 


100 


2.0x3.9 


104 


110 


88/11 


KP2M 


GoldCam+#240 


TI5 


4760- 


-5879 


1.52 


105 


2.0x3.9 


72 


111 


88/11 


MMTB 


BigBlue+300B 


Reticon 


3890- 


-7500 


1.34 


135 


2.5 


20 


112 


89/04 


MG24 


MarkIIIb+600V 


RCA 


5066- 


-6534 


2.91 


130 


1.7x10% 


31 


113 


89/06 


MG24 


MarkIIIb+600V 


TI4849 


5126- 


-6393 


2.18 


150 


2 4x10% 


23 


114 


89/06 


MMTB 


BigBlue+300B 


Reticon 


3890- 


-7500 


1.34 


135 


2.5 


12 


115 


89/08 


WHT4 


ISIS-Blue+R600B 


CCD-IPCS 


4330- 


-4970 


0.45 


95 


2.0x3.9 


8 


116^ 


89/10 


MMTR 


RedChannel+600B 


TI 


4750- 


-5950 


1.50 


125 


1.5x10% 


7 


117 


89/10 


MG24 


MarkIIIb+600V 


Thompson 


5031- 


-6300 


2.21 


130 


1 7x10% 


61 


118 


89/11 


MG24 


MarkIIIb+600V 


RCA 


5018- 


-6499 


2.89 


170 


1.7x 10% 


14 


119 


90/05 


MMTR 


RedChannel+600B 


TI 


4750- 


-5950 


1.50 


125 


1.5x10% 


17 


120 


90/10 


IN25 


IDS+235mm+R632V 


GEC6 


4806- 


-5606 


1.46 


100 


1.9x7.2 


10 


121 


91/05 


IN25 


IDS+235mm+R632V 


GEC3 


4806- 


-5603 


1.46 


100 


1.9x7.2 


87 


122 


91/10 


MG24 


MarkIIIb+600V 


Thompson 


5018- 


-6278 


2.20 


125 


1.7x10% 


43 


123 


91/11 


IN25 


IDS+235mm+R632V 


GEC3 


4806- 


-5603 


1.46 


100 


1.9x7.2 


29 


124 


92/01 


MG24 


MarkIIIb+600V 


Thompson 


5038- 


-6267 


2.15 


125 


1.7x10% 


35 


125 


92/06 


MG24 


MarkIIIb+600V 


Loral 


4358- 


-7033 


1.31 


125 


1.7x10% 


57 


126 


92/06 


CA22 


B&C spec+#7 


TEK6 


4800- 


-6150 


1.40 


100 


5.0x4.2 


39 


127 s 


92/09 


CT4M 


Argus+KPGL#3 


Reticon II 


3877- 


-6493 


2.19 


145 


1.9 


199 


128 


93/05 


MG24 


MarkIIIb+600V 


Loral 


5090- 


-7050 


1.40 


105 


1.7x10% 


2d 


129 


93/06 


MG24 


MarkIIIb+600V 


TEK 


4358- 


-6717 


2.31 


135 


1.7x10% 


3 


130^' 


93/06 


SS23 


DBS-blue+600B 


PCA 


5015- 


-5555 


0.80 


85 


2.0x10% 





131 s 


93/10 


ES36 


MEFOS+B&C+#26 


TEK512CB 


4850- 


-5468 


1.22 


105 


2.6 


128 


132 


93/10 


SS23 


DBS-blue+600B 


Loral 


4820- 


-5910 


1.10 


80 


2.0x10% 


14 


133 s 


93/10 


CT4M 


Argus+KPGL#3 


Reticon II 


3879- 


-6485 


2.19 


115 


1.9 


193 



a Telescopes: MG24=MDM 2.4m; KP4M/KP2M=KPNO 4m/2m; WHT4=William Herschel 4m (LPO); IN25=Isaac Newton 2.5m (LPO); 
MMTB/MMTR=MMT (blue/red); CA22=Calar Alto 2.2m; CT4M=Cerro Tololo 4m; SS23=Siding Spring 2.3m; ES36=ESO 3.6m. 
6 Spectral dispersion in A/pixel. 

Instrumental resolution (<r, not FWHM) in kms" 1 , as determined from the cross-correlation analysis calibration curves (see §3.1). 



d The aperture over which the galaxy spectrum was extracted: diameter for circular apertures and fibres, widthxlength for rectangular 
slits (10% means the spectrum was extracted out to the point where the luminosity had fallen to 10% of its peak value). 
e The number of spectra taken in the run. 
* These runs produced no useful data. 
9 These runs used fibre spectrographs. 



galaxy axes. The nominal goal in all cases was to obtain 
around 500 photons/ A at Mgfe, corresponding to a S/N per 
100 kms -1 resolution element of about 30. In fact our spec- 
tra have a median of 370 photons/ A at Mgb, corresponding 
to a S/N per 100 kms" 1 of 26 (see §|1|). 

In each run several G8 to K5 giant stars with known he- 
liocentric velocities were observed. These 'velocity standard 
stars' are used as spectral templates for determining red- 
shifts and velocity dispersions. In observing these standards 
care was taken to ensure that the illumination across the 
slit was uniform, in order both to remove redshift zeropoint 
errors and to mimic the illumination produced by a galaxy, 
thereby minimising systematic errors in velocity dispersion 
estimates. This was achieved in various ways: by defocussing 
the telescope slightly, by moving the star back and forth 
across the slit several times, or by trailing it up and down 
the slit. Such procedures were not necessary for standards 
obtained with fibre spectrographs, as internal reflections in 



the fibres ensure even illumination of the spectrograph for 
all sources. Very high S/N (typically > 10,000 photons/A) 
were obtained in order that the stellar templates did not 
contribute to the noise in the subsequent analysis. 

The normal calibration exposures were also obtained: 
bias frames, flatfields (using continuum lamps internal to 
the spectrographs or illuminating the dome) and spectra 
of wavelength calibration lamps before and/or after each 
galaxy or star exposure. In general we did not make use 
of spectrophotometric standards as fluxed spectra were not 
necessary and we wished to minimise overheads as much as 
possible. 

The calibration procedures were slightly different for 
the three large datasets taken using fibre-fed spectrographs 
at CTIO (runs 127 and 133) and ESO (run 131). Because 
of the need to calibrate the relative throughput of the fibres 
in order to perform sky subtraction, fibre observations al- 
ways included several twilight sky flatfield exposures. Each 
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velocity standard star was observed through several fibres 
by moving the fibres sequentially to accept the starlight. 



2.3 Reductions 

The reductions of both the longslit and fibre observations 
followed standard procedures as implemented in the IRAfQ 
MIDAS and Starlink Figaro software packages. We briefly 
summarise the main steps in the reduction of our longslit 
and fibre data below; further details can be found in Baggley 
(1996). 

The first stage of the reductions, common to all obser- 
vations, was to remove the CCD bias using a series of bias 
frames taken at the start or end of the night. These frames 
were median-filtered and the result, scaled to the mean level 
of the CCD overscan strip, was subtracted from each frame 
in order to remove both the spatial structure in the bias 
pedestal and temporal variations in its overall level. We also 
took long dark exposures to check for dark current, but in no 
case did it prove significant. Subsequent reductions differed 
somewhat for longslit and fibre observations. 

For longslit data, the next step was the removal of pixel- 
to-pixel sensitivity variations in the CCD by dividing by a 
sensitivity map. This map was produced by median-filtering 
the flatfield exposures (of an internal calibration lamp or 
dome lamp) and dividing this by a smoothed version of it- 
self (achieved by direct smoothing or 2D surface fitting) in 
order to remove illumination variations in the 'flat' field. If 
necessary (because of a long exposure time or a susceptible 
CCD), cosmic ray events were identified and interpolated 
over in the two-dimensional image using either algorithmic 
or manual methods (or both). 

The transformation between wavelength and pixel posi- 
tion in longslit data was mapped using the emission lines in 
the comparison lamp spectra. The typical precision achieved 
in wavelength calibration, as indicated by the residuals of 
the fit to the calibration line positions, was < 0.1 pixel, cor- 
responding to 0.1-0.3A or 5-15 kms -1 , depending on the 
spectrograph setup (see Table ^). The spectra were then 
rebinned into equal intervals of log A so that each pixel 
corresponded to a fixed velocity interval, Av = cAz = 
c(10 AlogA -l), chosen to preserve the full velocity resolution 
of the data. 

The final steps in obtaining longslit spectra are sky- 
subtraction and extraction. The sky level was measured from 
two or more regions along the slit sufficiently far from the 
target object to be uncontaminated by its light. To account 
for variations in transmission along the slit, the sky under 
the object was interpolated using a low-order curve fitted 
to the slit illumination profile. A galaxy spectrum was then 
extracted by summing along the profile, usually over the 
range where the object's luminosity was greater than ~10% 
of its peak value, but sometimes over a fixed width in arcsec 
(see Table ^). Standard star spectra were simply summed 
over the range along the slit that they had been trailed or 
defocussed to cover. 



* IRAF is distributed by the National Optical Astronomy Ob- 
servatories which is operated by the Association of Universities 
for Research in Astronomy, Inc. under contract with the National 
Science Foundation. 



For the fibre runs the individual object and sky spec- 
tra were extracted first, using a centroid-following algorithm 
to map the position of the spectrum along the CCD. The 
extraction algorithm fitted the spatial profile of the fibre, 
in order to remove cosmic ray events and pixel defects, and 
then performed a weighted sum over this fit out to the points 
where the flux fell to ~5% of the peak value. Next, the 
dome-illumination flatfield spectra were median-combined 
and a sensitivity map for each fibre constructed by divid- 
ing each fibre's flatfield spectrum by the average over all 
fibres and normalising the mean of the result to unity. The 
pixel-to-pixel variations in the CCD response were then re- 
moved by dividing all other spectra from that fibre by this 
sensitivity map. Wavelength calibration was accomplished 
using the extracted comparison lamp spectra, giving similar 
precision to the longslit calibrations, and the object spec- 
tra were rebinned to a log A scale. Using the total counts 
through each fibre from the twilight sky flatfield to give 
the relative throughputs, the several sky spectra obtained 
in each fibre exposure were median-combined (after manu- 
ally removing 'sky' fibres which were inadvertently placed on 
faint objects). The resulting high-S/N sky spectrum, suit- 
ably normalised to each fibre's throughput, was then sub- 
tracted from each galaxy or standard star spectrum. 

The final step in the reductions for both longslit and 
fibre data was to manually clean all the one-dimensional 
spectra of remaining cosmic ray events or residual sky lines 
(usually only the 5577A line) by linearly interpolating over 
affected wavelengths. 

2.4 Spectrum Quality 

We have two methods for characterising the quality of our 
spectra. One is a classification of the spectra into 5 quality 
classes, based on our experience in reducing and analysing 
such data. Classes A and B indicate that both the redshift 
and the velocity dispersion are reliable (with class A giving 
smaller errors than class B); class C spectra have reliable 
redshifts and marginally reliable dispersions; class D spectra 
have marginally reliable redshifts but unreliable dispersions; 
class E spectra have neither redshifts nor dispersions. The 
second method is based on the S/N ratio per 100 km s 
bin, estimated approximately from the mean flux over the 
restframe wavelength rang e used to determine the redshifts 
and dispersions (see §3.1) under the assumption that the 
spectrum is shot-noise dominated. These two measures of 
spectral quality are complementary: the S/N estimate is 
objective but cannot take into account qualitative problems 
which are readily incorporated in the subjective classifica- 
tions. Figure ^ shows example spectra covering a range of 
quality classes and instrumental resolutions. 

Figure ^ shows the S/N distribution for the whole sam- 
ple and for each quality class individually, and gives the total 
number of objects, the fraction of the sample and the me- 
dian S/N in each class. For the whole sample, 39% of the 
spectra have S/N>30, 70% have S/N>20, and 96% have 
S/N>10. The two quality measures are clearly correlated, 
in the sense that better-classed spectra tend to have higher 
S/N. However there is also considerable overlap in the S/N 
range spanned by the different classes. This overlap has var- 
ious sources: (i) factors other than S/N which affect the 
quality of the redshift and dispersion estimates, notably the 
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T — I — I — | — I 1 1 1 1 — I — I 1 — I — | 1 — 

A548-2_B 286109362 
a,= 100 S/N = 38 Q=A 
cz=1 1537 a=275±19 



T 



T 



T 



T 1 1 1 1 — I 1 — I — | — I 1 1 1 — | — I 

A2589_G 687133154 
a,= 145 S/N = 31 Q=A 
cz=12500 tr=165±14 





I I I | I I I I | I I I I | I I I I | I I I 

A496_2 194123421 
a r =100 S/N = 23 Q=B 
cz=9708 a=179±19 



A1983.G 368119331 
C7.= 125 S/N = 43 Q=A 
cz=13707 a=220±10 




I I I | I I I I | I I I I | I I I I | I I I 
J11_A 302105217 
ct,= 125 S/N = 20 Q=B 
cz = 8870 (7 = 226±31 



I I I | I I I I | I I I I | I I I I | I I I 
J12_A 310102159 
a,= 145 S/N=20 Q = B 
cz=12120 ct=250±40 






I I I | I I I I | I I I I | I I I I | I I I 
A397.C 153101074 
a,= 145 S/N=16 Q = C 
cz = 8961 a=191±26. 



A2197.G 568121334 
^=100 S/N=17 Q=C 
cz=8246 ct= 1 71 ±22 




1 



J14_D 346104318 
a,= 125 S/N=13 Q=C 
cz = 9256 cj = 213±55 






5200 5300 5400 5500 
Wavelength (A) 



5200 5300 5400 5500 
Wavelength (A) 



5200 5300 5400 5500 
Wavelength (A) 



Figure 1. Example spectra covering a range of quality classes and instrumental resolutions: the top, middle and bottom rows are spectra 
with quality classes A, B and C respectively; the left, central and right columns are spectra with resolutions 100, 125 and 145 km s -1 
respectively. The label for each spectrum gives the galaxy name, the GINRUNSEQ spectrum identifier, the instrumental resolution, the 
S/N and quality class of the spectrum, the redshift, the dispersion and its estimated error. Note that the panels show relative flux and 
have a false zero for viewing convenience. 



available restframe spectral range (which depends on both 
the spectrograph setup and the redshift of the target) and 
whether the object has emission lines; (ii) errors in estimat- 
ing the S/N (e.g. due to sky subtraction errors, the neglect of 
the sky contribution in computing the S/N for fainter galax- 
ies, or uncertainties in the CCD gain (affecting the conver- 
sion from counts to photons); (iii) subjective uncertainties 
in the quality classification, particularly in determining the 
reliability of dispersion estimates (i.e. between classes B and 
C). Both ways of determining spectral quality are therefore 
needed in order to estimate the reliability and precision of 
the spectroscopic parameters we measure. 

3 ANALYSIS 

3.1 Redshifts and Dispersions 

We derived redshifts and velocity dispersions from our spec- 
tra using the fxcor utility in IRAF, which is based on the 
cross-correlation method of Tonry & Davis (1979). We pre- 
ferred this straightforward and robust method to more elab- 
orate techniques since it is well-suited to the relatively mod- 
est S/N of our spectra. We used a two-step procedure, ob- 
taining an initial estimate of the redshift using the whole 



available spectrum and then using a fixed restframe wave- 
length range for the final estimates of redshift and veloc- 
ity dispersion. The procedure was applied in a completely 
uniform manner to all the spectra in our sample as far as 
differences in wavelength range and resolution would allow. 

The first step in the cross-correlation analysis is to fit 
and subtract the continuum of each spectrum in order to 
avoid the numerical difficulties associated with a dominant 
low-frequency spike in the Fourier transform. In the first 
pass through fxcor the continuum shape was fitted with a 
cubic spline with the number of segments along the spec- 
trum chosen so that each segment corresponded to about 
8000 km s -1 . Each iteration of the fit excluded points more 
than 1.5cr below or 3<r above the previous fit. In this way 
we achieved a good continuum fit without following broad 
spectral features. We then apodised 10% of the spectrum at 
each end with a cosine bell before padding the spectrum to 
2048 pixels with zeros. 

This continuum-subtracted, apodised spectrum was 
then Fourier transformed and a standard 'ramp' fil- 
ter applied. This filter is described by 4 wavenumbers 
(fci, /C2, fe, fci), rising linearly from to 1 between k\ and 
/c2 and then falling linearly from 1 to between hi, and ki- 
ln the first pass these wavenumbers were chosen to be fci=4- 
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Figure 2. The distribution of S/N with quality class. For each 
class the panels give the total number of spectra, the percentage 
of the whole sample and the median S/N. 

8 and fc2=9-12 (tailored to remove residual power from the 
continuum without affecting broad spectral features), and 
k3=N pix /3 and k4=N P i X /2 (N p i x is the number of pixels in 
original spectrum before it is padded to 2048 pixels; these 
choices attenuate high-frequency noise and eliminate power 
beyond the Nyquist limit at N P i X /2). The same procedures 
were also applied to the spectrum of the stellar velocity stan- 
dard to be used as a template. The cross-correlation of the 
galaxy and stellar template was then computed, and the 
top 90% of the highest cross-correlation peak fitted with a 
Gaussian in order to obtain a redshift estimate. 

This procedure was repeated for every template from 
that run, and the redshifts corrected to the heliocentric 
frame. Offsets in the velocity zeropoint between templates, 
measured as the mean difference in the redshifts measured 
with different templates for all the galaxies in the run, were 
typically found to be <;30kms -1 . These were brought into 
relative agreement within each run by choosing the best- 
observed K0 template as defining the fiducial velocity zero- 
point. Applying these offsets brought the galaxy redshifts 
estimated from different templates into agreement to within 
< 3 km s -1 . (The removal of run-to-run velocity offsets is de- 
scribed below.) The mean over all templates then gave the 
initial redshift estimate for the galaxy. 

This initial redshift was then used to determine the 
wavelength range corresponding to the restframe range 

=4770A to A ri'i a x 

=5770A. This range was chosen for use 
in the second pass through f xcor because: (i) it contains the 
MgI6 5174A band, H/3 486lA and the Fel 5207A and 5269A 



lines, but excludes the Nal D line at 5892A, which gives 
larger velocity dispersions than the lines in the region of Mgb 
(Faber & Jackson 1976); (ii) for redshifts up to our sample 
limit of cz=15000kms -1 this restframe wavelength range is 
included in the great majority of our spectra. The input for 
the second pass was thus the available spectrum within the 
range corresponding to restframe 4770-5770A. All but two 
of our runs cover the restframe out to at least 5330A for 
cz=15000 kms -1 ; the exceptions are run 115 (which is not 
used for measuring dispersions) and run 131 (which reaches 
restframe 5207A). 

In the second pass through fxcor we employed only 
minimal continuum subtraction based on a 1- or 2-segment 
cubic spline fit, preferring the better control over contin- 
uum suppression afforded by more stringent filtering at low 
wavenumbers. After considerable experimentation and sim- 
ulation, we found that the best filter for recovering velocity 
dispersions was a ramp with the same k% and k$ values as in 
the first pass, but with k2=0.Ql{N max — N m in), where N min 
and N max are the pixels corresponding to X m in and X max , 
and fci=0.75fc2- Again, the top 90% of the highest cross- 
correlation peak was fitted with a Gaussian. The position 
of this peak, corrected for the motion of the template star 
and the heliocentric motion of the earth relative to both the 
template and the galaxy, gave the final redshift estimate. 

The galaxy's velocity dispersion, a g , is in principle re- 
lated to the dispersion of the Gaussian fitted to the cross- 
correlation peak, o x , by a x = a'g + 2of (where Oi is the 
instrumental resolution; Tonry & Davis 1979). In practice 
this relationship needs to be calibrated empirically because 
of the imperfect match between the spectra of a broadened 
stellar template and a galaxy and the effects of the filter ap- 
plied to both spectra. The calibration relation between a x 
and (j g for a typical case is shown in Figure ^ (see caption 
for more details). We estimate the instrumental resolution 
for a given run from the mean value of the calibration curve 
intercepts for all the templates in the run (er^ ~ a x /\/2 when 
(jg=0); these are the values listed in Table 

The values of heliocentric radial velocity and veloc- 
ity dispersion were determined in this second pass through 
fxcor for each galaxy spectrum using all the templates in 
the same run. The final step is then to combine the redshift 
and dispersion estimates from each template, as summarised 
below. 

For the redshifts the steps involved were as follows: 

(i) Cases where the ratio of cross-correlation function peak 
height to noise (the R parameter defined by Tonry & Davis 
1979) was less than 2 were rejected, as were cases that dif- 
fered from the median by more than a few hundred kms -1 . 

(ii) The mean offset between the redshifts from a fiducial K0 
template and each other template was used to shift all the 
redshifts from the other template to the velocity zeropoint of 
the fiducial. These offsets were typically < 50 kms -1 . (iii) A 
mean redshift for each galaxy was then computed from all 
the unrejected cases using 2-pass 2a clipping, (iv) Any tem- 
plate which gave consistently discrepant results was rejected 
and the entire procedure repeated. The scatter in the red- 
shift estimates from different templates after this procedure 
was typically a few kms -1 . 

A very similar procedure was followed in combining ve- 
locity dispersions except that a scale factor rather than an 
offset was applied between templates: (i) Cases with R<4 
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Figure 3. A typical calibration curve showing the relation be- 
tween the width of the cross-correlation peak, a x , and the true 
velocity dispersion of the galaxy, cr g . The crosses are the individ- 
ual calibrations obtained by broadening each of the other tem- 
plates in the run and cross-correlating with the template being 
calibrated. The solid curve is the calibration curve used, a series of 
linear segments joining the median value of a x at each calibrated 
value of <7 3 . The dashed curve is the theoretical relation when no 
filtering is applied, = a* + 2c?, where ai is the instrumental 
resolution, in this case 145 km s" 1 . Note that the calibration curve 
flattens for a g < Cj, indicating that the true dispersion becomes 
increasingly difficult to recover as it drops below the instrumental 
resolution. 



were rejected, (ii) The mean ratio between the dispersions 
from a fiducial KO template and each other template was 
used to scale all the dispersions from the other template to 
the dispersion scale of the fiducial. These dispersion scales 
differed by less than 5% for 90% of the templates, (iii) A 
mean dispersion for each galaxy was then computed from 
all the unrejected cases using 2-pass 2a clipping, (iv) Any 
template with a scale differing by more than 10% from the 
mean was rejected as being a poor match to the programme 
galaxies and the entire procedure was then repeated. (Note 
that no significant correlation was found between scale fac- 
tor and spectral type over the range G8 to K5 spanned by 
our templates.) The scatter in the dispersion estimates from 
different templates after this procedure was typically 3-4%. 

Two corrections need to be applied to the velocity dis- 
persions before they are fully calibrated: (i) an aperture cor- 
rection to account for different effective apertures sampling 
different parts of the galaxy velocity dispersion profile, and 
(ii) a run correction to remove systematic scale errors be- 
tween different observing setups. The latter type of correc- 
tion is also applied to the redshifts to give them a com mon 
zerop oint . These two corrections are discussed below at 
and §3.5 respectively. 



3.4 



3.2 Linestrength Indices 

Once redshifts and velocity dispersions were determined, 
linestrength indices could also be measured using the pre- 
scription given by Gonzalez (1993). This is a refinement of 
the original 'Lick' system in which a standard set of bands 
was defined for measuring linestrength indices for 11 features 
in the spectra of spheroidal systems (Burstein et al. 1984). 
Gonzalez (1993), Worthey (1993) and Worthey et al. (1994) 
describe how this system has been updated and expanded 
to a set of 21 indices. Here we measure both the Mgb and 
Mg2 indices. 

The feature bandpass for Mg6 index is 5160. 1-5192. 6A, 
encompassing the Mg I triplet with components at 5166. 6A, 
5172. OA and 5183. 2A. The continuum on either side of 
the absorption feature is defined in bands covering 5142.6- 
5161.4A and 5191.4-5206.4A. Mgfe is an atomic index, and 
so is defined as the equivalent width of the feature in 
An gstroms, 



Mg6 = 



(1) 



where the integral is over the feature bandpass, S(X) is the 
object spectrum and C(A) is the linear pseudo-continuum 
defined by interpolating between two continuum estimates, 
taken at the midpoints of the blue and red continuum bands 
to be the mean values of the observed spectrum in those 
bands. 

Closely related to Mgb is the Mg2 index, for which 
the feature bandpass is 5154. 1-5196. 6A and the continuum 
bands are 4895.1-4957.6A and 5301.l-5366.lA. This index 
measures both the Mg I atomic absorption and the broader 
MgH molecular absorption feature. Mg2 is a molecular in- 
dex, and so is denned as the mean ratio of flux to local 
continuum in magnitudes, 



Mg 2 = -2.51o gl 



J S(X)/C(X)dX 



AX 



(2) 



where the integral is over the Mg2 feature bandpass, 
AA=42.5A is the width of that bandpass, and the pseudo- 
continuum is interpolated from the Mg2 continuum bands. 

In fact we will often find it convenient to express the 
Mg6 index in magnitudes rather than as an equivalent width. 
By analogy with the Mg2 index, we therefore define Mg6' to 
be 



Mg6' = -2.51og 10 1 



Mgb \ 
AA ) 



(3) 



where in this case AA=32.5A, the width of the Mgfe feature 
bandpass. 

In passing it should be noted that a different defini- 
tion of linestrength indices has sometimes been used (e.g. 
Worthey 1994, equations 4 and 5) in which the integral of 
the ratio of the object spectrum and the continuum in equa- 
tions |l] and ^ is replaced by the ratio of the integrals. This al- 
ternative definition has merits (such as simplifying the error 
properties of measured indices), but it is not mathematically 
equivalent to the standard definition. In practice, however, 
the two definitions generally give linestrengths with negligi- 
bly different numerical values. 

It is usual in studies of this sort to employ the Mg2 in- 
dex as the main indicator of metallicity and star-formation 
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history. However we find it useful for operational reasons to 
also measure the Mgb index. One problem is that the limited 
wavelength coverage of the spectra from some runs means 
that in a number of cases we cannot measure the Mg2 index 
(requiring as it does a wider wavelength range) although we 
can measure the Mgb index. We obtain Mgb for 676 objects 
(with 299 having repeat measurements) and Mg2 for 582 ob- 
jects (with 206 having repeat measurements) . Another prob- 
lem with Mg2 is that the widely-separated continuum bands 
make it more susceptible than Mgb to variations in the non- 
linear continuum shape of our unfluxed spectra, which result 
from using a variety of different instruments and observing 
galaxies over a wide range in redshift. We therefore present 
measurements of both Mgb and Mg2 : the former because it 
is better-determined and available for more sample galaxies, 
the latter for comparison with previous work. As previously 
demonstrated (Gorgas et al. 1990, j0rgensen 1997) and con- 
firmed here, Mgb and Mg2 are strongly correlated, and so 
can to some extent be used interchangeably. 

Several corrections must be applied to obtain a 
linestrength measurement that is calibrated to the standard 
Lick system. The first correction allows for the fact that the 
measured linestrength depends on the instrumental resolu- 
tion. Since all our spectra were obtained at higher resolution 
than the spectra on which the Lick system was defined, we 
simply convolve our spectra with a Gaussian of dispersion 
[a\ick ~ cr ?) 1//2 m order to broaden our instrumental resolu- 
tion tTi (see Table]]]) to the Lick resolution of 200 kms" 1 . 

The second correction allows for the fact that the mea- 
sured linestrength depends on the galaxy's internal velocity 
dispersion — a galaxy with high enough velocity dispersion 
cr 9 will have features broadened to the point that they extend 
outside their index bandpasses, and so their linestrengths 
will be underestimated. Moreover, if an absorption feature is 
broadened into the neighbouring continuum bands then the 
estimated continuum will be depressed and the linestrength 
will be further reduced. The 'cr-correction' needed to cali- 
brate out this effect can be obtained either by measuring 
linestrength as a function of velocity broadening for a set 
of suitable stellar spectra (such as the templates obtained 
for measuring redshifts and dispersions) or by modelling the 
feature in question. 

Although most previous studies have adopted the for- 
mer approach, we prefer to use a model to calibrate our 
indices, since we observe a dependence of the Mgb profile 
shape on a that is not taken into account by simply broad- 
ening stellar templates. Our simple model assumes Mgb to 
be composed of three Gaussians centred on the three Mg I 
lines at A 6 =5166.6A, A C =5172.0A and A r =5183.2A with 
corresponding relative strengths varying linearly with dis- 
persion from 1.0:1.0:1.0 at o^lOOkms" 1 to 0.2:1.0:0.7 at 
a=300 km s _1 . This dependence on dispersion is empirically 
determined and approximate (the relative strengths of the 
individual lines are not tightly constrained), but it does sig- 
nificantly improve the profile fits compared to assuming any 
fixed set of relative weights. Such variation of the Mgb profile 
shape reflects changes, as a function of velocity dispersion, 
in the stellar population mix and relative abundances (par- 
ticularly of Mg, C, Fe and Cr), which each affect the profile 
in complex ways (Tripicco & Bell 1995). 

Using the estimated value of the index to normalise the 
model profile and the effective dispersion (a\ + o"ii C fc) 1//2 



to give the broadening, we can estimate both the pro- 
file flux which is broadened out of the feature bandpass 
and the resulting depression of the continuum. Correct- 
ing for both these effects gives an improved estimate for 
the linestrength. Iterating leads rapidly to convergence and 
an accurate cr-correction for the Mgb and Mg2 indices. 
We find that the Mgb cr-correction is typically +4% at 
100 kms -1 and increases approximately linearly to +16% 
at 400 kms" 1 ; the Mg2 cr-correction is typically 0.000 mag 
up to 200 kms" 1 and increases approximately linearly to 
0.004 mag at 400 kms" 1 . 

Note that the usual method of determining the cr- 
correction by broadening standard stars ignores the depen- 
dence of profile shape on changes in the stellar population 
mix as a function of luminosity or velocity dispersion. Our 
tests indicate that by doing so, the usual method tends to 
overestimate Mgb for galaxies with large dispersions: by 2% 
at 200 kms" 1 , 6% at 300 kms" 1 and 14% at 400 kms" 1 . The 
two methods give essentially identical results for Mg2, since 
it has much smaller cr-corrections due to its wider feature 
bandpass and well-separated continuum bands. 

The other corrections that need to be applied to the 
linestrength estimates are: (i) an aperture correction to 
account for different effec tive apertures sampling different 
parts of the galaxy (§3.4); (ii) a run correction to remove 



systematic scale errors between different observing setups 
(§3.5); and (iii) an overall calibration to the Lick system 
determined by comparisons with literature data (§4.3). 



3.3 Error Estimates 

Error estimates for our redshifts, velocity dispersions and 
linestrengths come from detailed Monte Carlo simulations 
of the measurement process for each observing run. By cali- 
brating the errors estimated from these simulations against 
the rms errors obtained from the repeat mea surements that 
are available for many of the objects (see §3.6), we can obtain 
precise and reliable error estimates for each measurement of 
every object in our sample. 

The procedure for estimating the uncertainties in our 
redshifts and velocity dispersions was as follows. For each 
stellar template in each observing run, we constructed a 
grid of simulated spectra with Gaussian broadenings of 100- 
300 kms" 1 in 20 kms -1 steps and continuum counts corre- 
sponding to S/N ratios of 10-90 in steps of 10. For each 
spectrum in this grid we generated 16 realisations assum- 
ing Poisson noise. These simulated spectra were then cross- 
correlated against all the other templates from the run in 
order to derive redshifts and velocity dispersions in the stan- 
dard manner. The simulations do not account for spectral 
mismatch between the galaxy spectra and the stellar tem- 
plates, but for well-chosen templates this effect is only sig- 
nificant at higher S/N than is typically found in our data. 

Figure ^ shows the random error in redshift and the 
systematic and random errors in dispersion as functions of 
input dispersion and S/N for four of the larger runs. The 
systematic errors in redshift are not shown as they are neg- 
ligibly small (~lkms~ ), although the simulations do not 
include possible zeropoint errors. The systematic errors in 
dispersion are generally small (a few percent or less) for 
S/N>20, but become rapidly larger at lower S/N. The ran- 
dom errors in redshift increase for lower S/N and higher dis- 
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Figure 4. Redshift and dispersion errors as functions of input dispersion and S/N (labelling the curves) from the simulations of four 
of the larger runs. The top panel shows the random error in the redshift and the centre and bottom panels show the systematic and 
random error in the dispersion. The vertical dotted line indicates the instrumental dispersion of each run. 



persion, while the random errors in dispersion increase for 
lower S/N but have a broad minimum at around twice the 
instrumental dispersion. These curves have the general form 
predicted for the random errors from the cross-correlation 
method (Tonry & Davis 1979, Colless 1987). 

Given the dispersion and S/N measured for a spectrum, 
we interpolated the error estimates from the simulation for 
that particular observing run to obtain the systematic and 
random errors in each measured quantity. We used the re- 
sults of these simulations to correct the systematic errors 
in the velocity dispersions and to estimate the uncertain- 
ties in individual measurements of redshift and dispersion. 
For quality class D measurements of redshifts, where the 
spectra are too poor to estimate a dispersion and hence a 
reliable redshift error, we take a conservative redshift error 
of 50kms _1 . 

The linestrength error estimates were obtained by gen- 
crating 50 Monte Carlo realizations of the object spectrum 
with Poisson noise appropriate to the spectrum's S/N level. 
The Mgfe and Mg2 linestrengths were then measured for each 
of these realizations and the error estimated as the rms error 
of these measurements about the observed value. The error 
estimate obtained in this fashion thus takes into account the 
noise level of the spectrum, but does not account for errors 
in the linestrength due to errors in the redshift and disper- 
sion estimates, nor for systematic run-to-run differences in 
the underlying continuum shape. 

The estimated errors in the spectroscopic parameters 



are compared with, and calibrated to, the rms errors derived 
from repeat observations in §p^q. 



3.4 Aperture Corrections 

The velocity dispersion measured for a galaxy is the 
luminosity-weighted velocity dispersion integrated over the 
region of the galaxy covered by the spectrograph aperture. 
It therefore depends on (i) the velocity dispersion profile; 
(ii) the luminosity profile; (iii) the distance of the galaxy; 

(iv) the size and shape of the spectrograph aperture; and 

(v) the seeing in which the observations were made. In or- 
der to intercompare dispersion measurements it is therefore 
necessary to convert them to a standard scale. The 'aperture 
correction' this requires has often been neglected because it 
depends in a complex manner on a variety of quantities some 
of which are poorly known. The neglect of such corrections 
may account in part for the difficulties often found in rec- 
onciling dispersion measurements from different sources. 

The aperture correction applied by Davies et al. (1987) 
was derived by measuring dispersions for a set of nearby 
galaxies through apertures of 4"x4"and 16"xl6". In this 
way they used their nearby galaxies to define the velocity 
dispersion profile and obtained a relation between the cor- 
rected value, a cor , and the observed one, a a bs- This turned 
out to be an approximately linear relation amounting to a 
5% correction over the distance range between Virgo and 
Coma. 

More recently j0rgensen et al. (1995) have derived an 
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Figure 5. The distribution of aperture corrections to a standard 
metric aperture, (a) The distribution of the ratio of the observed 
metric apertures to the standard metric aperture (corresponding 
to 1.7 arcsec at the redshift of Coma), (b)-(d) The aperture cor- 
rections to this standard metric aperture for the dispersion, Mgft 
and Mg2 measurements. Note that a cor is the aperture-corrected 
dispersion and <j obs is the raw observed dispersion; likewise for 
the linestrengths. 



aperture correction from kinematic models based on data 
in the literature. Published photometry and kinematics for 
51 galaxies were used to construct two-dimensional models 
of the surface brightness, velocity dispersion, and rotational 
velocity projected on the sky. They found that the position 
angle only gave rise to 0.5% variations in the derived disper- 
sions and could thus be ignored. They converted rectangular- 
apertures into an 'equivalent circular aperture' of radius r ap 
which the models predicted would give the same dispersion 
as the rectangular slit. They found that to an accuracy of 
4% one could take r ap = 1.025(xy/n) 1 ^ 2 , where x and y are 
the width and length of the slit. 

From their models they then calculated the correction 
factor from the observed dispersion to the dispersion in some 
standard aperture. For a standard metric aperture, they 
found this aperture correction to be well approximated by a 
power law of the form 



(Jobs 



CZ \ 
CZ J 



(4) 



where a bs and cr cor are the observed and corrected disper- 
sions, ro is a standard aperture radius, defined to be 1.7 arc- 
sec, and czo is a standard redshift, defined as the redshift 
of Coma. The standard metric aperture is thus 0.54 h~ kpc 
in radius. Alternatively, one can correct to a standard rela- 
tive aperture (defined to be R e /8) using the same power law 
relation, 



O cor 
0"of>s 



He 



(5) 



This power law approximates the true relation to within 1% 
over the observed range of effective apertures (compare the 
distribution metric aperture sizes in Figure |i|a with Figure 4c 
of J0rgensen et al.). 

We also apply an aperture correction to our 



linestrengths. j0rgensen et al. noted that the radial gradient 
in the Mg2 index is similar to the radial gradient in log a, 
and so applied the same aperture correction for Mg2 as for 
log a. We adopt this procedure for Mg2. For Mgb we convert 
to Mgb' (Equation ^) and, assuming that the radial profile 
of Mgb' is similar to that of Mg2 (and hence log a), we apply 
the log a aperture correction to Mgb' before converting back 
to Mgb. 

The distributions of corrections to the standard met- 
ric aperture for the dispersions and linestrengths are shown 
in Figures |B|b-d. These corrections are generally positive, as 
most objects in our sample are observed through larger effec- 
tive apertures and are further away than J0rgensen et al.'s 
standard aperture and redshift. The corrections to standard 
relative apertures are quite similar, although having slightly 
greater amplitude and range. We choose to adopt the cor- 
rection to a standard metric aperture in order to minimise 
the size and range of the corrections and to facilitate com- 
parisons with dispersions and linestrengths in the literature. 



3.5 Combining Different Runs 

In comparing the redshifts, dispersions and linestrengths ob- 
tained from different runs we found some significant system- 
atic offsets. The origin of these run-to-run offsets is not fully 
understood. For the redshifts, the use of different velocity 
standard stars as the fiducials in different runs clearly con- 
tributes some systematic errors. For the dispersions, the cali- 
bration procedure we use should in principle remove instru- 
mental systematics; in practice, scale differences are com- 
mon, as is shown by the range of scale factors needed to 
reconcile velocity dispersions from various sources in the 
compilation by McElroy (1995; see Table 2). 

We cannot directly calibrate the measurements from 
each run to the system defined by a chosen fiducial run, 
as there is no run with objects in common with all other 
runs to serve as the fiducial. Instead, we use the mean off- 
set, A, between the measurements from any particular run 
and all the other runs. To compute this offset we separately 
compute, for each galaxy i, the error-weighted mean value 
of the measurements obtained from the run in question, Xij, 
and from all other runs, yik'- 



<*<> = 



<*) = 



E* 



(6) 



Here j runs over the nn observations of galaxy i in the target 
run and k runs over the rii observations of galaxy i in all 
other runs; Sij and Sik are the estimated errors in Xij and 
yik- We then take the average over all galaxies, weighting by 
the number of comparison pairs, to arrive at an estimate for 
the offset of the target run: 



A = 



Y, l -m i n i {{x x ) - ( yi )) 



E, 



mini 



(7) 



Here i runs over the I galaxies in the sample. We can reject 
outliers at this point by excluding galaxies for which the 
difference (xi) — (yi) is larger than some cutoff: for cz, log a, 
Mgb' and Mg2 we required differences less than 300kms _1 , 
0.2 dex, 0.1 mag and 0.1 mag respectively. The uncertainty, 
e, in this estimate of the run offset is given by 
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Table 2. Calibration of observing runs to a common system. 
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a Run 109 is the fiducial run for Mg6', denned to have zero offset. 
b Runs 104 and 129 have no objects in common with other runs. 
c Run 118 has no Mg2 measurements in common with other runs. 
d Runs 116 and 130 have no usable data. 
e Run 131 has no Mg2 measurements. 



(jTm^) 2 

where 8{xi) and 8{yi) are the error-weighted uncertainties 
in (xi) and (j/j) given by 

s(^) 2 - (E, *ti 3 ) -1 , W = (E, s., 2 )- 1 . (9) 

We subtract the offset determined in this manner from each 
run and then iterate the whole procedure until there are no 
runs with residual offsets larger than 0.5e. As a final step, 
we place the entire dataset (now corrected to a common 
zeropoint) onto a fiducial system by subtracting from all 
runs the offset of the fiducial system. Note that the run 
corrections for dispersion and Mgb are determined in terms 
of offsets in log a and Mgfe'. 

In order to maximise the number of objects with 
multiple measurements, we included the dataset from the 
'Streaming Motions of Abell Clusters' project (SMAC: 
M.J.Hudson, priv.comm.; see also Smith et al. 1997) in this 
analysis. There is a considerable overlap between the SMAC 
and EFAR samples which significantly increases the number 



of comparison observations and reduces the uncertainties in 
the run offsets. We chose to use the 'Lick' system of Davies 
et al. (1987; included in the SMAC dataset) as our fiducial, 
in order to bring the 7 Samurai, EFAR and SMAC datasets 
onto a single common system. This is not possible with Mgb, 
which is not measured in most previous work or by SMAC. 
We therefore chose run 109 (the Kitt Peak 4m run of Novem- 
ber 1988) as the Mg6 fiducial because it had a large number 
of high-quality observations and the systematics of the slit 
spectrograph are believed to be well understood. 

We checked that this procedure gives relative run cor- 
rections consistent with those obtained by directly compar- 
ing runs in those cases where there are sufficient objects in 
common. We have also compared our method with a slightly 
different method used by the SMAC collaboration to de- 
termine the run corrections for their own data and found 
good agreement (M.J.Hudson, priv.comm.). We carried out 
Monte Carlo simulations of the whole procedure in order to 
check the uncertainties in the offsets computed according 
to Equation ^. We found that this equation in general pro- 
vides a good estimate of the uncertainties, although when 
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the number of comparisons is small or involve a small num- 
ber of other runs it can under-estimate the uncertainties by 
up to 30%. Our final estimates of the uncertainties are there- 
fore derived as the rms of the offsets from 100 Monte Carlo 
simulations. 

Table [| lists the offsets for each run computed according 
to the above procedure, their uncertainties based on Monte 
Carlo simulations, the number of individual measurements 
(N) and the number of comparison pairs (N c ). Note that 
to correct our observed measurements to the fiducial sys- 
tem we subtract the appropriate run offset in Table ^ from 
each individual measurement. Of the 31 spectroscopic runs 
with usable data, only runs 104 and 129 have no objects in 
common with other runs and hence no run corrections; run 
118 has no Mg2 measurements in common and so no run 
correction for Mg2. 

Weighting by the number of individual measurements 
in each run, the mean amplitude of the corrections and 
their uncertainties are 28±8kms _1 in cz, 0.023±0.015 dex 
in logcr, 0.008±0.006 mag in Mg6' and 0.015±0.006 mag 
in Mg2. The significance of the individual run corrections 
(in terms of the ratio of the amplitude of the offset to its 
uncertainty) varies; however over all runs the reduced \ 2 
is highly significant: 15.7, 4.0, 3.3 and 11.4 for the correc- 
tions to the redshifts, dispersions, Mgb and Mg2 respectively. 
Application of the run corrections reduces the median rms 
error amongst those objects with repeat measurements from 
18kms _1 to 14kms~ 1 in redshift, 6.3% to 5.6% in disper- 
sion, 4.9% to 4.4% in Mg6 and 0.012 mag to 0.009 mag in 
Mg2. We also checked to see whether applying the run cor- 
rections reduced the scatter in external comparisons bet ween 
our data and measurements in the literature (see 4.3). We 



3 and §4.2) to match the rms errors from objects with 



found that although the scatter is dominated by the com- 
bined random errors, the corrections did reduce the scatter 
slightly in all cases. 

As another test of the run corrections for Mgfe' and Mg2 
(and also, more weakly, for logcr), we compared the Mgb'-<r 
and Mg2-cr distributions for each run (after applying the run 
corrections) with the global Mg&'-u and Mg2-<r relations 
derived in Paper V. Using the x 2 goodness-of-fit statistic to 
account both for measurement errors in the dispersions and 
linestrengths and for the intrinsic scatter about the Mg-<r 
relations, we find that for Mgb'-a there were two runs (113 
and 132) with reduced \ 2 greater than 3, while for Mg2-cr 
there was one such run (122). In all three cases the removal 
of 1 or 2 obvious outliers decreased the reduced x' 2 to a 
non-significant level. 



3.6 Calibrating the Estimated Errors 

Obtaining precise error estimates is particularly important 
because we will make extensive use of them in applying 
maximum likelihood methods to deriving the Fundamen- 
tal Plane and relative cluster distances for our sample. Al- 
though we have estimated the measurement errors as care- 
fully as possible, simulating the noise in the observations 
and the measurement procedures, some sources of error are 
likely to remain unaccounted-for and we may be systemati- 
cally mis-estimating the errors. We therefore auto-calibrate 
our errors by scaling the estimated errors in the combined 
measurements (the internal error estimate, based on the in- 
dividual measurement errors derived from simulations; see 



repeat measurements (an external error estimate). 

Figure ^| shows the differential and cumulative distri- 
butions of the ratio of the rms error to the estimated er- 
ror for each galaxy with repeat measurements of redshift, 
dispersion, Mg& and Mg2. The smooth curves are the pre- 
dicted differential and cumulative distributions of this ratio 
assuming that the estimated errors are the true errors. The 
top panel shows the comparison using the estimated errors 
(including all the corrections discussed above). For the red- 
shifts and linestrengths, the estimated errors are generally 
under-estimates of the true errors, since the ratio of rms to 
estimated errors tends to be larger than predicted. For the 
dispersions the estimated errors are generally over-estimates 
of the true errors, since this ratio tends to be smaller than 
predicted. For all quantities the assumption that the esti- 
mated errors are consistent with the true errors is ruled out 
with high confidence by a Kolmorogorov-Smirnov (KS) test 
applied to the observed and predicted cumulative distribu- 
tions. These differences between the estimated errors from 
the simulations and the rms errors from repeat measure- 
ments reflect the approximate nature of the S/N estimates 
and systematic measurement errors not accounted for in the 
simulations. 

In order to bring our estimated errors into line with the 
rms errors from the repeat measurements, we found it neces- 
sary to add 15 km s~ in quadrature to the estimated redshift 
errors, to scale the dispersion and Mg& errors by factors of 
0.85 and 1.15 respectively, and to add 0.005 mag to the Mg2 
errors. These corrections were determined by maximising the 
agreement of the observed and predicted distributions of the 
ratio of rms to estimated errors under a KS test (excluding 
outliers with values of this ratio >3.5). The corrections are 
quite well determined: to within a couple of kms -1 for the 
redshift correction, a few percent for the dispersion and Mg& 
corrections, and 0.001 mag for the Mg2 correction. Applying 
these corrections and repeating the comparison of rms and 
estimated errors gives the lower panels of Figure ^, which 
shows the good agreement between the rms errors from re- 
peat measurements and the calibrated errors estimates for 
the redshifts, dispersions and Mg linestrengths. 

The need for such a correction to the redshift errors may 
be due in part to the residual zeropoint uncertainties in the 
redshifts and in part to a tendency for the simulations to 
under-estimate the errors for high S/N spectra. The origin 
of the over-estimation of the dispersion errors is uncertain, 
although it may result from slightly different prescriptions 
for estimating the S/N in the observations and the simula- 
tions. The under-estimation of the linestrength errors may 
be due to neglecting the effects of errors in the redshift and 
dispersion estimates and the different continuum shapes of 
spectra from different runs when measuring linestrengths. 



4 RESULTS 

4.1 Individual Measurements 

The previous two sections describe the observations and 
analysis of our spectroscopic data. Table ^ lists the observa- 
tional details for each spectrum and the fully-corrected mea- 
surements of redshift, dispersion, Mgfe and Mg2, together 



© 0000 RAS, MNRAS 000, 000-000 



Peculiar motions of early-type galaxies - II 13 



CD 



60 - 



40 



20 




60 



Z 40 
E 

=5 



20 - 












1234 50 1234 50 1234 50 12345 
Acz rms /Acz est Aa rms /Aa est AMgb rms /AMgb est AMg2 rms /AMg2 est 



Figure 6. Comparison of the estimated errors derived from simulations and the rms errors for galaxies with repeat measurements of 
redshift, dispersion, Mgb and Mg2. Each panel shows the differential and cumulative distributions of the ratio of rms error to estimated 
error. The stepped curves are the observed distributions, while the smooth curves are the predicted distributions. The upper panels 
show the comparisons using the original estimated errors; the lower panels show the comparisons after correcting the estimated errors 
as described in the text. 



with their calibrated error estimates. Note that these error 
estimates are the individual measurement errors, and must 
be combined in quadrature with the run correction uncer- 
tainties given in Table |^ to give the total error estimate. 
We list the measurement errors rather than the total errors 
because the total errors are not independent, being corre- 
lated for objects in the same run. The version of the table 
presented here is abridged; the full table will be available 
upon publication from NASA's Astrophysical Data Center 
(ADC) and from the Centre de Donnees astronomiques de 
Strasbourg (CDS). 



The entries in Table q are as follows; Column 1 gives 
GINRUNSEQ, a unique nine-digit identifier for each spec- 
trum, composed of the galaxy identification number (GIN) 
as given in the master list of EFAR sample galaxies (Table 3 
of Paper I), the run number (RUN) as given in Table [l], and 
a sequence number (SEQ) which uniquely specifies the ob- 
servation within the run; column 2 gives the galaxy name, as 
in the master list of Paper I; column 3 is the telescope code, 
as in Table [j]; column 4 is the UT date of the observation; 
columns 5 & 6 are the quality parameter (with an asterisk if 
the spectr um s hows emission features) and S/N of the spec- 
trum (see §2.4); columns 7 & 8 are the fully-corrected helio- 
centric redshift cz (in km s" 1 ) and its measurement error; 
columns 9 & 10 are the fully-corrected velocity dispersion a 
(in kms -1 ) and its measurement error; columns 11 & 12 are 



column 15 provides comments, the full meanings of which 
are described in the notes to the table. 

There are 1319 spectra in this table. Note that 81 ob- 
jects from our original sample do not have spectroscopic 
observations and do not appear in the table (see the list 
of missing GINs in the table notes). Three of these are the 
duplicate objects (GINs 55, 435, 476) and three are known 
stars (GINs 131, 133, 191). Most of the others are objects 
which our imaging showed are not early-type galaxies, al- 
though there are a few early-type galaxies for which we did 
not get a spectrum. There are 34 spectra which are unusable 
(<3=E) either because the spectrum is too poor (13 cases) 
or because the object was mis-identified (20 cases) or is a 
known star (1 case, GIN 123). Of the 1285 usable spectra 
(for 706 different galaxies), there are 637 spectra with Q=A, 
407 with Q=B, 161 with Q=C and 80 with Q=D. 

The cumulative distributions of the total estimated er- 
rors in the individual measurements (combining measure- 
ment errors and run correction uncertainties in quadrature) 
are shown in Figure for quality classes A, B and C, and 
for all three classes together. The error distributions can be 
quantitatively characterised by their 50% and 90% points, 
which are listed in Table ^j. The overall median error in a 
single redshift measurement is 22 kms -1 , the median rela- 
tive errors in single measurements of dispersion and Mgfc are 
10.5% and 8.2%, and the median error in a single measure- 
ment of Mg2 is 0.015 mag. 



the fully-corrected Mg& linestrength index and its measure- 
ment error (in A); columns 13 & 14 are the fully-corrected 
Mg2 linestrength index and its measurement error (in mag) ; 
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Table 3. Individual spectroscopic measurements (abridged) 
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This is an abridged version of this table; the full table will be available upon publication from NASA's Astrophysical Data Center 
(ADC) and from the Centre de Donnees astronomiques de Strasbourg (CDS). The columns of this table give: (1) observation identifier 
(GINRUNSEQ); (2) galaxy name; (3) telescope used; (4) date of observation; (5) quality parameter; (6) signal to noise ratio; (7— 
8) rcdshift and estimated error; (9—10) velocity dispersion and estimated error; (11-12) Mgfe lincstrcngth and estimated error; (13—14) Mg2 
lincstrcngth and estimated error; and (15) notes on each observation. In the notes, 'double' means the EFAR galaxy is double; 'star' 
means the EFAR object is a star not a galaxy; 'mis-ID' means the spectrum is for some galaxy other than the nominated EFAR object; 
'mis-ID*' means the spectrum is for a nearby star rather than the EFAR object; 'Mgb at sky' means the object is at a redshift which 
puts Mgfe on the 5577A sky line; '#=#' notes the duplicated pairs in the EFAR sample (see Paper I; only the first of the two GINs is 
used); emission line objects (with an asterisk on Q) have the emission features listed; 'H/3 abs' or 'H/3 abs, [OIII]' means the redshift is 
based on the H/3 absorption feature (and [OIII] if present), as the spectrum stops short of Mgb (no dispersion or Mgb index is given for 
these objects). The objects for which we have no spectrum have GINs: 7, 20, 29, 30, 34, 35, 55, 62, 64, 67, 82, 83, 91, 104, 121, 131, 133, 
134, 161, 181, 191, 214, 225, 228, 231, 234, 256, 265, 309, 327, 391, 405, 407, 417, 434, 435, 442, 450, 451, 452, 458, 463, 464, 465, 470, 
475, 477, 483, 484, 486, 494, 516, 520, 521, 522, 523, 526, 544, 551, 553, 567, 569, 570, 575, 576, 577, 587, 594, 597, 603, 605, 624, 625, 
644, 671, 727, 760, 793, 798, 801, 901. 
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Table 4. The distribution of estimated errors per measurement 
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Figure 7. The cumulative distributions of estimated errors for 
individual measurements of redshift, velocity dispersion and Mgft 
linestrength. The distributions for quality classes A, B and C are 
shown as full, long-dashed and short-dashed lines respectively; 
the overall distribution is the thick full line, (a) The distribu- 
tion of estimated errors in cz; (b) estimated relative errors in cr; 
(c) estimated relative errors in Mgft; (d) estimated errors in Mg2. 



4.2 Combining Measurements 

We use a weighting scheme to combine the individual mea- 
surements of each quantity to obtain a best estimate (and its 
uncertainty) for each galaxy in our sample. The weighting 
has three components: 

(i) Error weighting: For multiple measurements Xi hav- 
ing estimated total errors Ai (the measurement errors and 
run correction uncertainties added in quadrature), we weight 
the values inversely with their variances, i.e. by A~ 2 . 

(ii) Quality weighting: We apply a weighting Wq which 
quantifies our degree of belief (over and above the estimated 
errors) in measurements obtained from spectra wit h diff cr- 

for 



ent quality parameters. Following the discussion in §2.4 
spectra with Q=A,B,C,D,E we use Wq=1, 1,1, 0.5,0 in com- 
puting redshifts, Wq=1, 1,0.5,0,0 in computing dispersions 
and Wq = 1,1,0.5,0,0 in computing linestrengths. 



(iii) Run weighting: we also apply a run-weigh ting 
Wr=0 to exclude run 115, for reasons explained in §2.1 
all other runs are given Wr=1. 

The combined estimate X is thus computed from the 
individual measurements Xi as the weighted mean 



(10) 



where Wi = A^ 2 WQiWm- The uncertainty in this weighted 
mean is computed as 



(EiWi) 



-1/2 



(11) 



This procedure is used to obtain combined estimates of the 
redshift, dispersion and linestrengths for each galaxy. We es- 
timate the overall quality Q as the highest quality amongst 
the individual measurements and obtain a combined esti- 
mate of the S/N as 



S / N = (T,i( S / N )i W Qi W Ki) 1/2 



(12) 



using the same weightings as for the dispersions (except 
when the overall quality is Q=D, when these weightings are 
omitted) . 

Table ^ gives the combined estimates of the spectro- 
scopic parameters for each galaxy in the EFAR sample. The 
version of the table presented here is abridged; the full table 
will be available upon publication from NASA's Astrophys- 
ical Data Center (ADC) and from the Centre de Donnees 
astronomiques de Strasbourg (CDS). The table lists: galaxy 
identification number (GIN), galaxy name, cluster assign- 
ment number (CAN; see §™, and the number of spectra, 
redshifts, dispersions and Mgfe and Mg2 linestrengths ob- 
tained for this object; then, for each of redshift, dispersion, 
Mg& and Mg2: the combined estimate, its estimated total er- 
ror (A) and the weighted rms error from any repeat observa- 
tions (S); finally, the combined S/N estimate and the overall 
quality parameter (with an asterisk if the galaxy possesses 
emission lines) . Note that only objects with useful measure- 
ments are included; hence the lowest quality class present in 
this table is Q=D, and the 7 galaxies with only Q=E spec- 
tra (GINs 123, 284, 389, 448, 599, 637, 679) in Table | are 
omitted. 

The cumulative distributions of uncertainties in the 
combined results are shown in Figure |^, both for the en- 
tire dataset and for quality classes A, B and C separately. 
The error distributions can be quantitatively characterised 
by their 50% and 90% points, which are listed in Table ^. 
The overall median error in redshift is 20kms~ 1 , the median 
relative errors in dispersion and Mg6 are 9.1% and 7.2%, and 
the median error in Mg2 is 0.015 mag. For the whole sam- 
ple, and for quality classes A and B, the median errors in 
the combined measurements are smaller than the median 
errors in the individual measurements, as one expects. How- 
ever for dispersion, Mg6 and Mg2 the errors are larger for 
quality class C and at the 90th percentile; this results from 
assigning a quality weighting of 0.5 to Q=C when combining 
the individual measurements of these quantities. 

The uncertainties listed in Table ^ represent the best 
estimates of the total errors in the parameters for each 
galaxy. However it must be emphasised that they are not 
independent of each other, as the run correction errors are 
correlated across all measurements from a run. To prop- 
erly simulate the joint distribution of some parameter for 
the whole dataset, one must first generate realisations of 
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Table 5. Spectroscopic parameters for the EFAR galaxies (abridged) 
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This is an abridged version of this tabic; the full table will be available upon publication from NASA's Astrophysical Data Center (ADC) 
and from the Centre de Donnees astronomiques dc Strasbourg (CDS). The columns of this table give: (1) galaxy identification number 
(GIN); (2) galaxy name; (3) the cluster assignment number (CAN); (4) the number of spectra N s , redshifts N z , dispersions N a , Mgb 
linestrcngths Nf, and Mg2 lincstrcngths 7V2 obtained for this object; then the combined estimate, its estimated total error (A) and the 
weighted rms error from any repeat observations (<5) for each of (5—7) rcdshift, (8—10) dispersion, (11-13) Mg6 linestrength and (14— 
16) Mg2 linestrength; (17) the combined S/N estimate; and (18) the overall quality parameter (with an asterisk if the galaxy possesses 
emission lines). Only objects with useful measurements are included; hence the lowest quality class present in this tabic is Q=D, and the 
7 galaxies with only Q=E spectra (GINs 123, 284, 389, 448, 599, 637, 679) arc omitted. 
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Table 6. The distribution of errors per galaxy 
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Figure 8. The cumulative distributions of the total estimated 
errors in the combined measurements of redshift, velocity disper- 
sion, Mgfe and Mg2 for each galaxy. The distributions for quality 
classes A, B and C are shown as full, long-dashed and short- 
dashed lines respectively; the overall distribution is shown as the 
thick full line, (a) The distribution of combined errors in cz; 
(b) combined relative errors in u; (c) combined relative errors 
in Mgb; (d) combined errors in Mg2. 



the run correction errors (drawn from Gaussians with stan- 
dard deviations given by the uncertainties listed in Table ^|) 
and the individual measurement errors (drawn from Gaus- 
sians with standard deviations given by the uncertainties 
listed Table ||). For each individual measurement, one must 
add the realisation of its measurement error and the realisa- 
tion of the appropriate run correction error (the same for all 
measurements in a given run) to the measured value of the 
parameter. The resulting realisations of the individual mea- 
surements are finally combined using the recipe described 
above to yield a realisation of the value of the parameter for 
each galaxy in the dataset. 

The distributions of redshift, velocity dispersion, Mgb 
and Mg2 for the galaxies in the EFAR sample are displayed 
in Figure ^| The galaxies for which we measured velocity 
dispersions are only a subset of our sample of programme 
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Figure 9. The distributions of (a) redshift, (b) velocity disper- 
sion, (c) Mgft linestrength, and (d) Mg2 linestrength for the galax- 
ies in the EFAR sample. 
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Figure 10. The fraction of programme objects for which we mea- 
sured a velocity dispersion as a function of the logarithm of the 
selection diameter Dw (in arcsec). 



galaxies (629/743), and represent a refinement of the sam- 
ple selection criteria. Figure shows the fraction of pro- 
gramme galaxies with measured dispersions as a function of 
the galaxy diameter Dw on which the selection function of 
the programme galaxy sample is defined. There is a steady 
decline in the fraction of the sample for which usable dis- 
persions were measured, from 100% for the largest galaxies 
(with Dw > 40 arcsec) to about 75% for the smallest (with 
8 arcsec <J Dw <, 15 arcsec; there are only 3 programme 
galaxies with Dw < 8 arcsec). This additional selection 
effect must be allowed for when determining Fundamental 
Plane distances. 



4.3 Internal and External Comparisons 

One of the strengths of our spectroscopic sample is the high 
fraction of objects with repeat observations: there are 375 
galaxies with a single dispersion measurement, 160 with two 
measurements and 141 with three or more measurements. 
Figure [n] shows the cumulative distributions of rms errors 
in redshift, dispersion, Mgb and Mg2 obtained from these re- 
peat observations. The detailed internal comparisons made 
possible by these repeat measurements have been used to 
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Figure 11. The cumulative distributions of the rms errors from 
repeat measurements of redshift, velocity dispersion, Mgfe and 
Mg2. The distributions for spectral types A, B and C are shown as 
full, long-dashed and short-dashed lines respectively; the overall 
distribution is shown as the thick full line, (a) The distribution of 
rms errors in cz in kms" ; (b) relative rms errors in it; (c) relative 
rms errors in Mgb; (d) rms errors in Mg2. 



establish the run cor rections (§3.5) and to calibrate the es- 
timated errors (§3.6). The latter process ensured that the 



estimated errors were statistically consistent with the rms 
errors of the repeat measurements. 

We also make external comparisons of our measure- 
ments with the work of other authors. The EFAR redshifts 
are compared in Figure |l^ with redshifts given in the litera- 
ture by the 7 Samurai (Davies et al. 1987), Dressier & Sheet- 
man (1988), Beers et al. (1991), Malumuth et al. (1992), 
Zabludoff et al. (1993), Colless & Dunn (1996) and Lucey 
et al. (1997). Only 11 of the 256 comparisons give redshift 
differences greater than 300 kms" 1 : in 6 cases the EFAR 
redshift is confirmed either by repeat measurements or other 
published measurements; in the remaining 5 cases the iden- 
tification of the galaxy in question is uncertain in the lit- 
erature. For the 245 cases where the redshift difference is 
less than 300 kms , there is no significant velocity zero- 
point error and the rms scatter is 85 kms -1 . Since our repeat 
measurements show much smaller errors (90% are less than 
36 kms -1 ), most of this scatter must arise in the literature 
data, some of which were taken at lower resolution or S/N 
than our data. 
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Figure 12. Differences between EFAR redshifts and those from 
various sources in the literature. 



Figure |l3| compares the EFAR dispersions with pub- 
lished dispersions from the work of the 7 Samurai (Davies 
et al. 1987), Guzman (1993), j0rgensen et al. (1995) and 
Lucey et al. (1997), and the compilation of earlier measure- 
ments by Whitmore et al. (1985). Note that we do not com- 
pare to the more recent compilation by McElroy (1995), 
since its overlap with our sample is essentially just the sum 
of above sources. The mean differences, A = logcr_BF.4fl — 
log ant, and their standard errors are indicated on the fig- 
ure; none of these scale differences is larger than 6% and 
in fact all five comparisons are consistent with zero scale 
error at the 2a level or better. The rms scatter in these 
comparisons is significantly greater than the errors in our 
dispersion measurements, implying that in general the lit- 
erature measurements have larger errors and/or that there 
are unaccounted-for uncertainties in the comparison. 

We determine the zeropoint calibration of our 
linest reng th measurements with respect to the Lick system 
(see §3.2) by comparing our Mgb' and Mg2 linestrengths 
to measurements for the same galaxies given by Trager 
et al. (1998). We find that slightly different calibrations 
are needed for objects with different redshifts, the re- 
sult of slight variations in the non-linear continuum shape 
as the spectra are redshifted with respect to the instru- 
ment response and the sky background (see §3.2). Good 



agreement with Trager et al. is obtained if we use differ- 
ent zeropoints for galaxies with redshifts above and be- 
low cz=3000 kms -1 (although there are no objects in the 
comparison at cjz>10000 kms -1 ). Excluding a few out- 
liers, we find weighted mean differences between the EFAR 
and Trager et al. linestrengths of (AMgr/}=— 0.022 mag 
and (AMg 2 ) =-0.083 mag for cz<3000 kms" 1 , and 
(AMgfo')=-0.008 mag and (AMg 2 }=-0.028 mag for 
cz>3000 kms -1 . Subtracting these zeropoint corrections 
gives the final, fully-corrected, linestr eng th measurements 
as listed in Tables |5| & ^. Figures |lj & |l5f show the residual 
differences between the EFAR and Trager et al. linestrength 
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Figure 13. Comparisons of EFAR dispersions with those 
from various sources in the literature: (a) Davies et al. (1987), 
(b) Guzman (1993), (c) j0rgensen (1997), (d) Lucey et al. (1997), 
and (e) Whitmore et al. (1985). In each case the mean difference, 
A = (log &EFAR — logcr; it ), and its standard error are indicated, 
along with the rms scatter and the number of galaxies in the 
comparison. 
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Figure 14. The residual differences in Mgb' lincstrcngths from 
EFAR and Trager et al. (1998) after applying the zeropoint cor- 
rections discussed in the text: (a) the distribution of residuals; 
(b) the residuals as a function of Mgb'; (c) the residuals as a 
function of logo - ; (d) the residuals as a function of redshift. Out- 
liers excluded from the determination of the zeropoint corrections 
are shown as crosses. 



measurements after applying these zeropoint corrections. 
The rms scatter is 0.019 mag in Mgi/ for the 41 objects 
in common, and 0.023 mag in Mg2 for the 24 objects in 
common. There is no statistically significant trend with 
linestrength, velocity dispersion or redshift remaining in the 
residuals after these zeropoint corrections are applied. 

Figure ^ compares our calibrated Mg2 linestrengths to 
those obtained in A2199, A2634 and Coma by Lucey et al. 
(1997). The overall agreement for the 36 objects in common 
is very good, with a statistically non-significant zeropoint 
offset and an rms scatter of 0.029 mag, similar to that found 
in the comparison with Trager et al. 

The relation between the measured Mg6' and Mg2 
linestrengths for all the galaxies in the EFAR sample is 
shown in Figure [l7| We fit this relation using a maximum 
likelihood technique which accounts for both measurement 
errors and selection effects (Saglia et al. 1998, in prepara- 
tion; Paper VI). We find 

Mg 2 = 1.94Mg6' - 0.05 , (13) 

with a perpendicular rms residual of 0.019 mag (correspond- 
ing to an rms of 0.041 mag in Mg2, or 0.021 mag in Mg6'). 
The relation is the same if we fit ellipticals, E/SOs, cDs or 
spirals separately. This relation is similar to those derived by 
Burstein et al. (1984) and j0rgensen (1997). We can there- 
fore use Mg6' as a predictor of Mg2 (albeit with larger un- 
certainties) for those cases where Mg2 cannot be measured 
directly. 

Also shown in Figure [It] is the predicted relation be- 
tween Mgfo' and Mg2 as a function of age and metallicity 
given by Worthey (1994). His models correctly predict the 
slope of the relation, but are offset by —0.025 mag in Mgb'(or 
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Figure 15. The residual differences in Mg2 linestrengths from 
EFAR and Trager et al. (1998) after applying the zeropoint cor- 
rections discussed in the text: (a) the distribution of residuals; 
(b) the residuals as a function of Mg2 ; (c) the residuals as a func- 
tion of log a\ (d) the residuals as a function of redshift. Outliers 
excluded from the determination of the zeropoint corrections are 
shown as crosses. 
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Figure 16. Comparisons of Mg2 linestrengths obtained by EFAR 
and Lucey et al. (1997). The mean difference, A = Mg 2 (EFAR) 
— Mg2 (Lucey), and its standard error are indicated, along with 
the rms scatter and the number of galaxies in the comparison. 
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Figure 17. The relation between Mgfe' and Mg2 and its max- 
imum likelihood fit. Ellipticals are marked by circles, E/SOs by 
squares, cDs by asterisks and spirals by triangles. Typical esti- 
mated errors are shown in the top left corner. The relation be- 
tween Mgb' and Mg2 as a function of age and metallicity, as pre- 
dicted by Worthey (1994), is shown as the grid lying parallel to, 
but offset from, the data. 



by +0.05 mag in Mg'2), indicating a difference in the model's 
zeropoint calibration for one or both indices. 



5 CLUSTER ASSIGNMENTS 

The correct assignment of galaxies to clusters (or groups) is 
crucial to obtaining reliable redshifts and distances for the 
EFAR cluster sample. We also need to increase the precision 
of the cluster redshifts in order to minimise uncertainties in 
the clusters' peculiar velocities. To achieve these goals we 
merged the EFAR redshifts with redshifts for all galaxies in 
ZCAT (Huchra et al., 1992; version of 1997 May 29) which 
lie within 3h _1 Mpc (2 Abell radii) of each nominal EFAR 
cluster centre (see Table 1 of Paper I). We then examined 
the redshift distributions of the combined sample in order 
to distinguish groups, clusters and field galaxies along the 
line of sight to a nominal EFAR 'cluster'. We also considered 
the distribution of galaxies on the sky before assigning the 
EFAR galaxies to specific groupings. 

The results of this process are shown in Figure |ll| which 
shows the redshift distributions of galaxies within 3h _1 Mpc 
around each of the nominal EFAR clusters (labelled by their 
cluster ID number, CID; see Paper I) and the adopted group- 
ings in redshift space. Note that CID=81 (A2593-S) does 
not appear since it was merged with CID=80 (A2593-N) — 
see below. Each EFAR galaxy was assigned to one of these 
groupings and given a cluster assignment number (CAN), 
listed in Table H. The main grouping along the line of sight 
has a CAN which is simply the original two-digit CID; other 
groupings have CANs with a distinguishing third leading 
digit. The groupings (which we will hereafter call clusters 
regardless of their size) are labelled by their CANs in Fig- 
ure ^ which also shows the boundaries of each cluster in 
redshift space. The last two digits of each galaxy's CAN is its 
CID, apart from 41 galaxies which were reassigned to other 
neighbouring clusters: two galaxies in CID=33 were reas- 
signed to CAN=34 (GINs 254, 255); two galaxies in CID=34 
were reassigned to CAN=33 (GINs 263, 264); five galaxies 
in CID=35 were reassigned to CAN=36 (GINs 270, 274, 
275, 281, 282); fourteen galaxies in CID=36 were reassigned 
to CAN=35 (GINS 285-292, 295-297, 299-301), one galaxy 
in CID=47 was reassigned to CAN=50 (GIN 406); three 
galaxies in CID=59 and two in CID=61 were reassigned to 
CAN=53 (GINs 514, 517, 527, 536, 537); five galaxies with 
CID=69 were reassigned to CAN=70 (GINs 617, 618, 619, 
622, 623); and all seven galaxies with CID=81 were reas- 
signed to CAN=80 (GINs 709-715). 

Table (?] lists, for each CAN, the number of EFAR galax- 
ies, the number of EFAR+ZCAT galaxies, and the mean red- 
shift, its standard error (taken to the error in the redshift 
for clusters with only one member) and the velocity disper- 
sion. These quantities are computed both from the EFAR 
sample and from the EFAR+ZCAT sample. In many of the 
clusters the EFAR sample is greatly supplemented by the 
ZCAT galaxies, leading to much-improved estimates of the 
mean cluster redshift: using EFAR galaxies only the median 
uncertainty in the mean cluster redshift (for clusters with 
more than one member) is 177kms _1 ; with EFAR+ZCAT 
galaxies the median uncertainty is reduced to 133kms _1 . 
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Figure 18. The redshift distributions of galaxies within 3h _1 Mpc of each nominal EFAR cluster using the EFAR and ZCAT data. 
Each distribution is labelled at top right by the nominal cluster ID number (CID). The solid histogram shows the distribution of EFAR 
galaxies; the open histogram shows the extra ZCAT galaxies. The groupings adopted have boundaries in redshift marked by dotted lines 
and are labelled by their cluster assignment number (CAN). Clusters without numbers and boundaries contain no EFAR galaxies. 
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Table 7. Cluster mean redshifts and velocity dispersions 

EFAR EFAR+ZCAT EFAR EFAR+ZCAT . 
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6 CONCLUSIONS 

We have described the observations, reductions, and anal- 
ysis of 1319 spectra of 714 early-type galaxies studied as 
part of the EFAR project. We have obtained redshifts for 
706 galaxies, velocity dispersions and Mg6 linestrengths for 
676 galaxies, and Mg2 linestrengths for 582 galaxies. Al- 
though obtained in 33 observing runs spanning seven years 
and 10 different telescopes, we have applied uniform proce- 
dures to derive the spectroscopic parameters and brought all 
the measurements of each parameter onto a standard sys- 
tem which we ensure is internally consistent through com- 
parisons of the large numbers of repeat measurements, and 
externally consistent through comparisons with published 
data. We have performed detailed simulations to estimate 
measurement errors and calibrated these error estimates us- 
ing the repeat observations. 

The fully-corrected measurements of each parame- 
ter from the individual spectra are given in Table ^| 
the final parameters for 706 galaxies, computed as the 
appropriately-weighted means of the individual measure- 
ments, are listed in Table [| The median estimated errors 
in the combined measurements (including measurement er- 
rors and run correction uncertainties) are Acz—20 kms -1 , 
A(t/<t=9.1% (i.e. Alogcr=0.040 dex), AMgfc/Mgfc=7.2% 
(i.e. AMgfc'=0.013 mag) and AMg 2 =0.015 mag. Compar- 
isons with redshifts and dispersions from the literature show 
no systematic errors. The linestrengths required only small 
zeropoint corrections to bring them onto the Lick system. 

We have assigned galaxies to physical clusters (as op- 
posed to apparent projected clusters) by examining the line- 
of-sight velocity distributions based on EFAR and ZCAT 
redshifts, together with the projected distributions on the 
sky. We derive mean redshifts for these physical clusters, 
which will be used in estimating distances and peculiar ve- 
locities, and also velocity dispersions, which will be used to 
test for trends in the galaxy population with cluster mass 
or local environment. 

The results presented here comprise the largest single 
set of velocity dispersions and linestrengths for early-type 
galaxies published to date. These data will be used in com- 
bination with the sample selection criteria of Wegner et al. 
(1996, Paper I) and the photometric data of Saglia et al. 
(1997, Paper III) to analyse the properties and peculiar mo- 
tions of early-type galaxies in the two distant regions studied 
by the EFAR project. 
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spectrograph with grating #7 (60 A/mm) was used in com- 
bination with the TEK#6 CCD. 

ESO 3.6m: The MEFOS fibre feed and the Boiler & 
Chivens spectrograph were used. MEFOS has 58 2.6 arcsec 
diameter fibres (29 for targets and 29 for sky) positioned 
within a 1 degree diameter field at prime focus. The detector 
was a Tektronix TK512CB CCD (ESO#32). 

CTIO 4m: The ARGUS 24-object fibre spectrograph 
was used. ARGUS has a 50 arcmin field at the f/2.8 prime 
focus. Each of the 24 arms holds two 1.9 arcsec diameter 
fibres which lie 36 arcsec apart on the sky; one arm is po- 
sitioned on the target and the other on sky. The fibres feed 
a thermally and mechanically isolated bench spectrograph 
with a 510mm focal length Schmidt blue collimator and a 
229mm focal length Schmidt camera. A Reticon If 1200x400 
CCD detector was used with grating KPGL #3. 



APPENDIX A: OBSERVING DETAILS 

This appendix gives further details of the instrumental con- 
figurations used on different telescopes. 

MDM 2.4m: The Mark Ilia spectrograph, with a 
1.87 arcsec wide slit, was used for all runs up to the end 
of 1988; from 1989 this was replaced by the Mark Illb 
spectrograph, which is identical except that a 1.68 arcsec 
slit was used (except for run 113, when the slit width was 
2.36 arcsec). For runs 101-103 a 600 lines/mm grism blazed 
at 4600A was used; for all subsequent runs, a 600 lines/mm 
grism blazed at 5700 A was employed. The slit was usually 
oriented N-S. Two-pixel binning perpendicular to the dis- 
persion direction was employed to lower the readout noise. 

KPNO 4m: The RC spectrograph and grating KPC-17B 
(527 lines/mm) were used with the UV Fast Camera and the 
TI2 CCD. 

KPNO 2.1m: The Gold spectrograph/camera and grat- 
ing #240 (500 lines/mm) were used with the TI5 CCD. 

WHT 4.2m: The blue arm of the ISIS spectrograph 
was used with the CCD-IPCS imaging photon counting sys- 
tem. Most objects were observed using the R600B grating 
(600 lines/mm), but one object (J26 A, GIN=648) was ob- 
served with the R300B grating (300 lines/mm). 

INT 2.5m: The Intermediate Dispersion Spectrograph 
(IDS) and R632V grating (632 lines/mm) were used with 
the 235mm camera for all runs. 

SSO 2.3m: Both runs used the blue arm of the Double 
Beam Spectrograph (DBS) with a 600 lines/mm grating. 
Run 130 used the Photon Counting Array (PCA) , while run 
132 used a Loral CCD. 

MMT Blue: The 'Big Blue' spectrograph was employed 
with a 300 lines/mm grating (blazed at 4800 Ain first order) 
and the Reticon detector. The MMT image stacker gave two 
2.5 arcsec circular apertures separated by 36 arcsec. 

MMT Red: The MMT Red Channel was used with a 
600 lines/mm grating (blazed at 4800 A) and the 800x800 TI 
CCD binned by two pixels perpendicular to the dispersion. 
The slit was 1.5 x 180 arcsec, but heavily vignetted in the 
outer 30 arcsec in one direction. 

Calar Alto 2.2m: The Cassegrain Boiler & Chivens slit 
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